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We discuss algebraic/numeric methods to compute one-loop corrections for multiparticle/jet production cross 
sections. By using efficient reduction algorithms a compact expression for the ggg"/"/ amplitude is obtained. 
Further a numerical approach for 6-point 1-loop diagrams is presented. 



1. INTRODUCTION 

The theoretical description of multi-particle 
production at the one-loop level is a very chal- 
lenging task, as the complexity of the Feynman 
diagrammatic approach grows exponentially with 
the number of external partons. No Standard 
Model process which has generic 2^4 kinemat- 
ics is computed at the one-loop level although this 
is highly relevant for many Higgs boson search 
channels at the LHC, like giuon fusion and weak 
boson fusion, where additional jets have to be 
tagged to improve the signal to background ratio. 
For signal reactions like PP H + 0,1,2 jets, 
with H jj,WW* ,T^T^ which are available 
at one-loop level, many backgrounds remain to be 
calculated. As an example for needed calculations 
consider PP bbbb + X, PP ^ + 2 jets -I- X 
or PP — > ZZ -f 77 -f AT, which require the eval- 
uation of hexagon graphs like the ones given in 
Fig. [2 

Te computation of the related amplitudes relies 
on efficient methods for the evaluation of the cor- 
responding Feynman graphs. In the next section 
we shortly review our reduction formalism. As 
an example for the efficiency of our methods we 
discuss the 5-point 1-loop amplitude gg 773 
in Section 3. It seems to be feasible to apply the 
presented techniques also for 6-point processes, 
as long as the internal masses of the problem can 
be neglected. In {1<2) we have shown that in the 
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Figure 1. Typical Hexagon graph for the 1-loop 
amplitudes gg — > bbbb and gg ZZjj. 



massless Yukawa model our formalism leads to 
compact expressions. Going to the massive case 
leads generally to much more involved expressions 
and in that case numerical methods seem to be 
preferable. An approach for the numerical evalu- 
ation of 6-point Feynman diagrams is outlined in 
Section 4. 

2. REDUCTION FORMALISM 

In the Feynman diagrammatic approach any 
one-loop amplitude can be represented as a linear 
combination of factors which contain the group 
theoretical information and tensor one-loop inte- 
grals: 

where 
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The propagator momenta are defined by qj = fc — 
Tj = k — pi . . . — pj. To separate the Lorentz 
structure from the integrals it is useful to express 
the tensor integrals in terms of scalar integrals 
with nontrivial numerators. 
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[i?/2] is the smaller nearest integer to R/2. The 
bracket with the Lorentz indices as superscripts 
stands for the sum over all distinguishable dis- 
tributions of the Lorentz indices to the metric 
tensors and external momenta. The separation 
of the kinematical information allows to sort the 
amplitude into gauge invariant subsets. Basic in- 
gredient of the given formula is the Feynman pa- 
rameter integral defined in _D = n + 2™ dimen- 
sions: 

/^(ji,...,jfl)-(-irr(7v-^/2) 

(1) 

In 12] we have derived a reduction formula for such 
parameter integrals for general R in arbitrary 
dimensions D. It is based on differentiation by 
parts in parameter space. The derived formula 
maps rank R A^-point integrals in D dimensions 
to rank R — 1 A^-point integrals and higher di- 
mensional integrals with lower rank. As the lat- 
ter are IR finite, a separation of IR divergent and 
IR finite terms can be obtained in this way which 
defines an approach for a semi-numeric method. 
After extraction of all UV/IR poles the remaining 
integrals can be treated numerically. If one wants 
to proceed analytically one has to iterate the for- 
mula. In this way one can show that arbitrary 



iV-point Feynman integrals can be expressed in 
terms of n = 4 — 2e dimensional bubble and trian- 
gle functions and ?i + 2 dimensional boxes. More 
details on reduction formalisms can be found in 
IB4. 



3. THE LOOP AMPLITUDE gg 
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To give an example for our algebraic approach 
we have considered the 5-point 1-loop amplitude 
gg 5 . This amplitude is indirectly known 

from the 1-loop 5-giuon amplitude |S] by turning 
gluons into photons. 

We define all particles as incoming. 



7(pi, Ai) + 7(p2, A2) + 5(p3, A3, C3) 
+5(P4, A4, C4) + 5(p5, A5, C5) -> 



(2) 



In hadronic collisions this amplitude is relevant 
for the production of photon pairs in associa- 
tion with a jet and as such a contribution of the 
background to the Higgs boson search channel 
H 77 -I- jet. For a phenomenological analysis 
see (7I8| . The colour structure of this amplitude 
can be written as 

T^^^^'^^-f^/ggg ^ 0] = :i9^/'=3C4C5_4Ai A3A4A5 

y^AiA2A3A4A5 helicity dependent linear combi- 
nations of scalar integrals and a constant term 
which is a remnant of two-point functions with 
coefficients of order (D — 4). Six independent he- 
licity components exist: +++++,++++ -, 

- ++++, \-++, +++ — As the 

amplitude is finite one expects that all 3-point 
functions which carry spurious infrared poles can- 
cel. The function basis of the problem is thus 
reduced to 2-point functions 

r(i + e)r(i-e)2 i-s,,)-^ 

2^^^'^- r(2-2e) e ' 

4-point functions in 6 dimensions written as [2] 
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and constant terms. From unitarity one expects 
that the -I— I— I— I— I- , -I— I— I— I- - , - amplitudes 
should be polynomial. The other helicity ampli- 
tudes will also contain non-polynomial functions 
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like logarithms and dilogarithms contained in 
and Fi. 

To give an example for a compact helicity am- 
plitude we show here the result for A 
only. The remaining ones which have also com- 
pact representations can be found in The re- 
sult is expressed in terms of field strength tensors 
J^!^'^ = Pj^.'j ~ P'j^j where are the polarization 
vectors of the gluons and photons. 

We split the result of A ^^'^ into three pieces 
with indices F,B,1, which belong to the part pro- 
portional to 6-dimcnsional boxes -Fi, a part con- 
taining bubble graphs , and a constant term, 
respectively. 

A-+++ = A-p-+++ + As-+++ + A^-+++ 
We find 



Cp- +++ Pi ■ T+ ■ P3 -(3^4) Fi(si3, Si4, 525) 
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The indicated permutations have to be applied to 
the indices of the field strength tensors, momenta 
and Mandelstam variables. The coefficients are 
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In the given expressions the S2 ® S3 symmetry 
under exchange of the two photons and the three 
gluons is manifest after taking into account the 
omitted colour factor. 

The result indicates that with our approach 
indeed a compact representation of complicated 
loop amplitudes can be obtained. The application 
of our approach to relevant 6-point amplitudes is 
presently under study. 

4. NUMERICAL APPROACH 

Due to the complexity of the analytic approach 
if massive particles are present, a numerical ap- 
proach seems to be more appropriate to tackle 
different types of one-loop amplitudes in a uni- 
fied and efficient way. 

Recently a great activity in that direction with 
many new ideas can be observed |9lll)lllll^ . 

4.1. Reduction to basic building blocks 

As basic building blocks for an amplitude in 
our numeric approach, we choose scalar 2-point 
functions and 3-point functions /J and n + 2 
dimensional box functions 14 with nontrivial 
numerators. The latter are infrared finite. Possi- 
ble UV singularities are only contained in the 2- 
point functions and their subtraction is straight- 
forward. The (soft and coUinear) IR singularities 
are, as a result of the reduction, only contained in 
2-point functions and 3-point functions with one 
or two light-like legs. In this form, they are easy 
to isolate and to subtract from the amplitude. 
After reduction and separation of the divergent 
parts, we are left with finite integrals ^3 (ji, J2, J3) 
and ji, 3a), with nontrivial numera- 

tors. As numerical stability problems are entirely 
from the denominators we discuss only the case 
of scalar integrals with trivial numerators here. 
Systematic methods for the combination of the 
IR divergences from the virtual corrections with 
their counterparts from the real emission contri- 
bution already exist (^Sl ^tnd references therein). 

In this section we focus on the evaluation of a 
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finite 6 point scalar integral. As a first step we 
reduce the hexagon integral to box and triangle 
functions which are the basic building blocks of 
the reduction. 

4.2. Parameter representation of basic 
building blocks 

To evaluate the box and triangle functions nu- 
merically, we first perform a sector decomposi- 
tion. 

1 = Q{xi > X2, ■ ■ ■ ,xn) + Q{X2 > Xi, . . . ,xn) 
+ . . . + Q{xn > xi,. . . ,xn^i) (3) 

for the integration over N parameters {N — 3 for 
the triangle, N = A for the box). The step func- 
tion is defined as 1 if the inequality of its ar- 
gument is fulfilled, and else. Now, we carry out 
one parameter integration explicitly. We show 
the explicit expressions only for the triangle inte- 
gral, the ones for the box are analogous and can 
be found in We obtain 

tD / 2 2 2\ 

^3 {si,S2,S3,mi,m2,m^) = 

oD I 2 2 2\ 

I aD I 2 2 2\ 



(si,S2,S3,Wi,m2,m3) 



with 



= / dtidt2 



A 



1 



{l + ti+ t2) Atl +Bt2 + C-iS 



(4) 



B = {rrii + 7712 ~ ^2)^1 + '772 + '^3 ~ S3 
C = m\t\ + (777^ 4- TTig — Si)ti + m\ 
R = B'^ - 4:AC + i6 
T^2A{l + ti)-B 

The discussion is also valid in the case of vanish- 
ing masses or invariants, as long as the functions 
remain IR finite. Note that if infrared divergences 
are present the triangle integrals can typically be 
treated analytically. The (r7-|- 2)-dimensional box 
function are infrared finite for any physically rel- 
evant kinematics. 



4.3. Singularity structure 

Starting from Q one integration is performed 
explicitly. In order to analyse the singularity 
structure of the integrands, we then separate 
imaginary and real part. One obtains 



^Tti (Sl, 52, S3, 7771,7772, 7773) 
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log (\B ~ y/R\j - iTr&{B <^/R<2A + B 
T + ^/R 



^ ■flogf|2A + B + \/;R 

-log(|S + /R|) +77re(B < < 2yl + B)) } 

Three regions which lead to an imaginary part 
can be distinguished: 

Region I: A + B + C -2A < B < 0, 

C > <^ (E < ±y/R <2A + B). 

Region II: A+B+C > 0, C < ^ {B < VR < 

2A + B) and not {B < -^R <2A + B). 

Region III: A + B + C <Q,C >Q ^ 

{B < -/R < 2A + B) and not [B < ^fR < 
2A + B). 

Region I is an overlap region where the imaginary 
part has two contributions. In regions II and III 
only one of the 0-functions contributes. Note 
that the box function /f has the same singu- 
larity structure [H]. As /f and /f are the 
basic building blocks, this analysis of the singu- 
larity structure is done once and for all. Knowing 
the critical region of integration it is possible to 
map out the singularities by adequate parameter 
transformations. 
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4.4. Numerical integration 

To demonstrate the practicality of our method 
to evaluate multi-leg integrals, we show in Fig. |21 
a scan of the 2mj = 350 GeV threshold of the 
4-dimensional scalar hexagon function for a real- 
istic kinematical configuration. For details of the 
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Figure 2. Scan of the 2mt — 350 GeV threshold of 
the 4-dimensional scalar hexagon function which 
corresponds kinematically to the right Feynman 
diagram of Fig. 



integration methods see |14I15) . 

5. CONCLUSION 

To make reliable phenomenological studies for 
collider experiments operating at the TeV scale 1- 
loop calculations with many external particles are 
mandatory. In this talk I have outlined recent de- 
velopments concerning the analytic and numeric 
evaluation of 1-loop Feynman diagrams. Using 
reduction methods a compact result for the 3- 



gluon 2-photon amplitude was presented. Con- 
cerning numerical methods we have developed 
an approach to successfully integrate hexagon 
functions numerically. Merging and applying 
these techniques to more challenging situations 
is presently under study. 
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